Observational constraints on phantom-like braneworld cosmologies 
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We investigate a simple braneworld model in which the universe contains only cold dark matter and 
a cosmological constant, but the effective dark energy is phantom-like because of extra-dimensional 
gravity effects. Modified gravity screens the cosmological constant A, allowing for a larger A. In 
practice, observations do not favour any significant screening. We use supernova data, the cosmic 
microwave background shift parameter, and the baryon oscillation peak in the galaxy distribution 
to constrain the model. We find the mean value of S7 m with 68% confidence limits, and an upper 
limit on Qa at the 68% confidence level. The best-fit model is very close to a standard LCDM 
model, but the LCDM model provides a better fit since it has one less parameter. 



I. INTRODUCTION 

There is increasingly strong evidence for the accel- 
eration of the late-time universe, from observations of 
supernovae, cosmic microwave background anisotropics 
and the large-scale structure. This late-time accelera- 
tion poses a major theoretical challenge in cosmology 0. 
Within the framework of general relativity, the accelera- 
tion originates from a dark energy field with effectively 
negative pressure: w =p/p < — |. (It is also conceivable, 
though it remains controversial, that acceleration maybe 
due to nonlinear backreaction or averaging effects 0].) 
The simplest model, with cosmological constant or vac- 
uum energy (w = —1) as dark energy, i.e., the LCDM 
model, provides a very good fit to the data 0, but the 
unnaturally small and fine-tuned value of A cannot be 
explained within current particle physics. Quintessence 
models, i.e. scalar fields with w > — 1, allow for richer 
dynamical behaviour, but do not improve the fit to the 
data, and also do not lessen the severe theoretical prob- 
lems faced by LCDM. "Phantom" models, with w < — 1, 
violate the null energy condition p+p > 0, and as a result, 
the total phantom energy density grows with expansion. 
Phantom scalar fields share the theoretical problems of 
quintessence, but in addition, they have negative kinetic 
energy and lead to an unstable quantum vacuum 

Current observations are compatible with w < — 1. 
The WMAP 3-year data, in combination with large-scale 
structure and SN data, allows for w < —1 in a general 
relativistic model: for constant w |3j, 

(1) 



w = -1.06] 



The simplest way to produce w < — 1 is a phantom scalar 
field in general relativity. But the price is instability. 
Thus it is interesting to investigate other models in which 
w < —1, but without negative kinetic energy or associ- 
ated instabilities. This can happen if some other effect 
mimics phantom dynamics, i.e., leads to an effective w e s 
such that w c g < — 1, but without the presence of a patho- 
logical phantom field. 

For example, Boisseau et al. [j| showed that this can 
happen in scalar-tensor theories. A number of other pos- 



sible mechanisms for mimicking phantom behaviour has 
recently been investigated , including nonminimal cou- 
pling to gravity, interactions between dark energy and the 
matter sector, modified Lagrangians, braneworld mod- 
els and quantum effects. Many of these mechanisms re- 
quire complicated additional features and fine-tunings - 
it is not easy to construct simple and natural phantom- 
like behaviour without a phantom field. Among the less 
complicated models is that of Csaki et al. 7], in which 
photon-to-axion conversion mimics super-acceleration (in 
the presence of non- phantom dark energy). 

Phantom behaviour may also occur in braneworld 
models with an infrared modification of general relativ- 
ity, as pointed out by Sahni and Shtanov Q. A special 
case of the Sahni-Shtanov models was further investi- 
gated by Lue and Starkman : 9]. In these models, the 
4-dimensional brane universe contains only matter and 
a cosmological constant A, but a 5-dimensional gravita- 
tional effect leads to phantom behaviour. The most im- 
portant implication of this modified gravity effect is that 
A is effectively screened, so that in principle a higher 
value of A is allowed by the observations than in the 
LCDM case. These models are the simplest in the Shani- 
Shtanov class, with only one parameter more than the 
LCDM model, and they are the simplest braneworld 
models with phantom behaviour as far as we are aware. 
In this paper, we extend the investigation of Lue and 
Starkman by testing the models against observations, to 
see how the data constrains the model parameters, and 
to compare the goodness of fit to that of LCDM. As ex- 
plained below, this test can only be partially carried out, 
since the CMB anisotropies and matter power spectrum 
for the model have not yet been computed, and this re- 
mains a formidable open problem. However we can apply 
geometric observational tests that are based on the back- 
ground dynamics. 



II. THE PHANTOM-LIKE BRANEWORLD 

The Dvali-Gabadadze-Porrati (DGP) braneworld 
model, generalized to cosmology by Deffayet 0] i is a seuC_ 
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accelerating model without any form of dark energy. In 
fact, the self-accelerating models are not the only form of 
DGP model. There are two separate branches, DGP(±), 
depending on how the 4D brane universe is embedded in 
the 5D spacetime. 

• The (+) branch is the self-accelerating model, 
which has recently been tested against super- 
nova (SN) and baryon oscillation (BO) observa- 
tions [llj, and then in addition against the CMB 
shift data [l^- The flat DGP(+) model is outside 
the 2(7 contour for the joint SN and BO constraints, 
but this is somewhat misleading, since the CMB ob- 
servations have not been applied. When the CMB 
shift constraint is included, the flat DGP(+) model 
is within the 2ct contour for the joint constraints, 
but the best-fit LCDM model gives a significantly 
better fit to the data [l2| . 

• The DGP(— ) model is very different. It does not 
self- accelerate, but requires dark energy on the 
brane. The simplest model has a cosmological con- 
stant, and we will call this the LDGP model, fol- 
lowing Ref. 9j. The LDGP model is a special case 
of the Sahni-Shtanov models with zero brane ten- 
sion. It experiences 5D gravitational modifications 
to its dynamics, which effectively screen the cos- 
mological constant. At late times, as gravity leaks 
off the 4D brane, the dynamics deviates from gen- 
eral relativity. The transition from 4D to 5D be- 
haviour is governed by a crossover scale r c , as in the 
(+) branch. (The LDGP model has recently been 
generalized by replacing the cosmological constant 
with a quintessence field [lflj-) 

The energy conservation equation for LDGP remains 
the same as in general relativity, but the Friedman equa- 
tion is modified: 



p + 3H(p + p) = 0, 



H 



H 8nG 



A 
3~ 



(2) 
(3) 



[The self-accelerating (+) branch has —H/r c instead of 
+H/r c .] For the CDM case, with p = 0, these equations 
imply 



H 



-AirGp 



y/1 + 32nGr^p/3 + 4r^A/3 



(4) 



Equation J2J shows that at early times, the general rela- 
tivistic Friedman equation is recovered: 



H > rZ 



8ttG A 
— p+ 3 



(5) 



By contrast, at late times, the H/r c term is important 
and the Friedman equation is nonstandard. 



Gravity leakage at late times screens the cosmological 
constant, leading to an effective dark energy 



where 



H 2 = —-(p + p eS ). 



By Eq. |@$, H < 0, so that 

Peff > . 



(6) 



(7) 



(8) 



This is the basis for the phantom-like behaviour in the 
model. We can define the effective dark energy equation 
of state w Q g = Peg/peg via 

p cS + 3H(l + w cS )p cS = 0, (9) 

or equivalently, 

H = -ATvG[p+{l+w eS )p eS \. (10) 

Thus p e fi and w e s, as defined in Eqs. (j5|)- (jTU)l . are ef- 
fective quantities that give a standard general relativis- 
tic interpretation of LDGP expansion history, i.e., they 
describe the equivalent general relativity model. This 
equivalent model is a phantom model: by Eqs. 10 and 
@, we have 



1 + w cS 



H 



8irGr c Hp c s 



(11) 



Since H < by Eq. (|10|l . we have effective phantom 
behaviour, 



w cS < -1 (/Ccff > 0) 



(12) 



The effective picture breaks down when p e g = 0, i.e., 
when H = = r c A/3, which always happens at some 
redshift z* in the history of the universe. At this time, 
u> c ff — ► —00. At higher redshifts, 1 + w e s is positive, de- 
creasing from +00, as illustrated in Fig.^ This is a spe- 
cial case of the behaviour found by Sahni and Shtanov p| • 
Nothing unphysical occurs at z* , it is only the description 
that breaks down. Note that it is possible that < 0, 
i.e., in the future, so that Hq > H*. But in this case, 
u> e ff at z = is positive, so that there is no phantom be- 
haviour for nonnegative redshifts. In fact, Eq. 1221) below 
shows that this only occurs for the unphysical parameter 
values, O m > 1. 

Although w c ff — > — 5 for large redshifts, CDM is be- 
coming strongly dominant, so that the geometrical con- 
straints from observations can be comfortably satisfied 
- as we show in the next section. As discussed further 
in the concluding section, it remains an open question 
whether structure formation constraints can also be sat- 
isfied. 



(a) 



(b) 





FIG. 1: (a) and (b): 1 + w e g against 1 + z. (c): 1 + w to t against 1 + z, with the same parameters as (a). 

(a) shows a typical model with w e « < — 1 at z = 0, with parameters (f2 mi flA) = (0.3,0.8). (b) shows a model with unrealistic 
parameters, (f2 m ,fiA) = (1.1, 1.4), that has w e s > —1 at z = 0. Note that w e s — » —1 as z — > —1 (a — > oo), and w e s — » — 5 as 
« — » 00. The redshift 2* where iu e fi goes through ±00 is clearly evident, with 2* > in (a) and 2* < in (b). Note that w to t is 
continuous at 2*. 



The effective phantom behaviour in LDGP has no asso- 
ciated instability, unlike phantom scalar fields in general 
relativity. Furthermore, in general relativistic phantom 
models, H eventually becomes positive, i.e., the universe 
eventually super-accelerates, which can lead to a "big 
rip" singularity |l4j . This also happens in a number of 
models that mimic phantom behaviour without phantom 
matter. By contrast, in LDGP, H is always negative, and 
there is no big rip singularity in LDGP. Equation l|10l) 
shows that the key issue is the sign of 1 + u> e ff as a — > oo. 
In phantom and most phantom-like models, 1 + w e g < 
as a — > 00, so that H becomes positive. In LDGP, we 
have 1 + w c ff — > 0~ and the universe is asymptotically de 



Sitter: 



Weft 



-1 , H -> H, 



dS 



1 

2^ 



4r?A 



1 - 1 



(13) 



Note that H^s is les s th an the corresponding LCDM 
value -HdS.LCDM = \/A/3. This is an aspect of the ef- 
fective screening of A by DGP(— ) gravity. 

In LDGP, the total equation of state parameter is al- 
ways greater than — 1, i.e., the phantom effects never 
dominate. Although the screened dark energy has phan- 
tom equation of state, w e g < — 1, the total w remains 
above —1, and there is no phantom-like acceleration of 
the universe, and thus no big rip. This is shown by 
Eq. I|23|) below, and illustrated in Fig. ^ 
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III. OBSERVATIONAL CONSTRAINTS ON 
LDGP 



In dimensionless form, the LDGP Friedman equa- 
tion J3J) implies 

e{z) = ^- =^n ni (i + z)* + n A + n ra -. 



where 



i2m 



SirGpo 
1 

4Hf! 



A 
3H§ 



Note that Eq. © implies 

Q m + Qa - 
so that in particular, 

fi m + Ov > 1 



(15) 
(16) 

(17) 
(18) 



This condition ensures that H% > 0. The flat LCDM 
model is the limiting case of equality, when Q rr = 0. 
The unphysical region f2 m + < 1 (corresponding to 
imaginary Hq) is shown in Fig. 

Equation (|17|l means that the flat LDGP model has 
the remarkable feature that it mimics a closed LCDM 
model in the (On>Ov) plane: 



Ok 
2 



'Or. <-> — ^ ■ (19) 
The WMAP 3-year data with SN data gives Q 

n K = -o.oii;^ , n A = 0.72 ± .04 . (20) 



If we take the mean value, Vt^ = —0.01, then Eq. H19|) 
gives the "equivalent" value r c = IOOHq 1 . This very 
large value of r c suggests that observations will require 
the LDGP model to be close to the LCDM model (r c = 
oo). We confirm this expectation below. 
By Eqs. lfTT|) and (|T4*|) . we have 



1 + Weff(z) 



/ a~Q m (l + z) 3 



E{z)] 



(21) 



This equation shows explicitly how we can mimic the ex- 
pansion history of LDGP by an effective phantom model 
in general relativity. However, there is no violation of the 
null energy condition, and no super-acceleration, since in 
LDGP the phantom dynamics is gravitational. Figure ^ 
illustrates this equation. 
At the current time, 



1 + w cfr (0) = - 



1)0 



(l-OnXOn + QA + l) 



(22) 



Provided that fl m < 1, we have phantom behaviour at 
the current time, i.e., Eq. (|12fl is satisfied. 



As noted above, the phantom effects never domi- 
nate, and the total equation of state parameter, w to t = 
Ptot/ptot = w c spcs/(p + Pes), is always greater than -1. 
By Eqs. and (gjj , 



Wtot(z) 



On(l + zf 



E(z) 



E(z)] ' 



(23) 



and it follows immediately that Wtot{z) > —1. Since E — > 
Q A + fl,, c — v/OT asz — > — 1 (i.e. , as a — * oo) , it follows 
that 1 + Wtot - ► + . At early times, i.e., z — > oo, we 
have w t ot ~~ > 0, which is the GR limit (we are neglecting 
radiation). This behaviour is illustrated in Fig. 2] 
From Eq. (@}, the dimensionless acceleration is 



a/a 

fl 



E 



2n rc + 2fl A - Q m (l + zf 



(24) 



2 V 'On(i + z) 3 + n A + O c 

so that the redshift when acceleration starts is given by 



2<M 

oj 



1/3 



>9l 




1/3 



(25) 



The LCDM result is recovered for f2 rc 
0, this gives the critical line f^A = 2f2 n 

On < -(1 + A ), 



0. When z acc 
- 1, so that 



(26) 



is the condition for models that are currently accelerat- 
ing. This differs from the corresponding condition for 
LCDM, i.e., fl m < 2tt A . 

The fundamental test of the background dynamics of 
a cosmological model is the SN magnitude-redshift test, 
based on the luminosity distance, 



d L = 



(1 + *) 



dz' 

W) 



(27) 



(Here we are restricting to the flat case.) The 68%, 95% 
and 99% confidence contours from fits to the Legacy 
(SNLS) data [H are shown in the LDGP parameter 
plane in Fig. 

Further independent tests are needed to check whether 
the LDGP model is consistent with the observed features 
of the universe. 

The CMB shift parameter 



S — V On-ffo 



dh{Zr) 
(1 + Zr) 



(28) 



encodes the relation between the angular diameter dis- 
tance to last scattering, the angular scale of the first 
acoustic peak, and the physical scale of the sound hori- 
zon. This parameter is effectively model-independent 
and provides a good test of the background dynamics, 
independent of the SN redshift test. We take z r = 1090. 
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FIG. 2: The 68%, 95% and 99% confidence contours for the SNLS data (left), the CMB shift parameter S (based on WMAP 
3-year data) (middle) and the baryon oscillation measure A (right). The shaded area is the unphysical region, Q m + Qa < 1, 
and the boundary is the flat LCDM model, Q rc — 0. 



n A 1 




FIG. 3: Joint constraints on LDGP models from the SNLS 
data, the BO measure A and the CMB shift parameter S. 
The solid dot is the best-fit model. The shaded area is the 
unphysical region, f2 m + Qa < 1, and the boundary is the flat 
LCDM model, n rc = 0. 



Wang and Mukherjee 0] have used the WMAP 3-year 
data to find that 



S= 1.70 ±0.03. 



(29) 



This value shows a significant improvement in the error 
over the value from the 1-year data 0], S = 1-72 ±0.06. 
The constraints from the CMB S parameter, based on 
the 3-year data, are shown in Fig. |2] 

The baryon acoustic oscillation peak recently detected 
in the SDSS luminous red galaxies (LRGs) can 
be used as another independent geometrical test of the 



LDGP. The correlation function for SDSS LRGs shows 
a peak at a scale ~ 100/i -1 Mpc, corresponding to the 
first acoustic peak at recombination (determined by the 
sound horizon). The observed scale effectively constrains 
the quantity |ig 



a = Jn 



H$dl( Zl ) 



-,1/3 



(30) 



where Z\ = 0.35 is the typical LRG redshift. Eisenstein 
et al. 0| give 



A = 0.469 ± 0.017 . 



(31) 



(We have suppressed a weak dependence of A on the 
spectral tilt.) Note that there is a level of uncertainty in 
the use of the BO measure A to test non-LCDM mod- 
els. This uncertainty will require substantial further ef- 
fort, both in the processing of the data fl9| and in 
the theoretical understanding of the LDGP matter power 
spectrum (see Ref. an d the discussion below). For 
the present, we use the constraint in Eq. (|31|l . and as- 
sume that the unresolved issues will lead to only small 
corrections. Constraints from the BO data are shown in 
Fig. El 

The joint constraints from SN, CMB shift and BO data 
are shown in Fig. [3J 



IV. CONCLUSIONS 

A naive comparison between LDGP models and closed 
LCDM models, as illustrated in Eqs. JTHJ) and J2DJ|, al- 
ready suggests that the data will favour LDGP models 
that are close to the flat LCDM model. The joint con- 
straint contours in Fig. [21 confirm this expectation. 

The best fit from the joint constraints is 



fin 



0.26, ft A = 0.74, w eS (0) 



-1.00. 



(32) 
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FIG. 4: The normalized marginalized likelihoods for fi m and J1a, subject to the joint constraints of data from SNLS supernova 
redshifts, CMB shift parameter S, and baryon acoustic oscillation measure A. 



with x 2 of 

X 2 = 114.6 , x 2 per degree of freedom = 0.996 . (33) 

Interestingly, the best-fit LDGP model, with current data 
and within the accuracy of our calculations, is the LCDM 
limit of LDGP. Since the LCDM model has one less pa- 
rameter, it provides a better fit as measured by the re- 
duced x 2 - 



X 



LCDM 



per degree of freedom = 0. 



(34) 



It is the CMB shift value from the 3-year WMAP data 
that is forcing the best-fit strongly towards the LCDM 
limit: if we use the shift value from the 1-year data [l7j . 
i.e., S = 1.72 ± 0.06, then the best-fit model from the 
joint constraints moves off the LCDM boundary: 

fi m = 0.27 , n A = 0.78 , Wcfi(O) = -1.01 (1-year S) , 

(35) 

with x 2 = 113.7, and a reduced \ 2 P er degree of freedom 
of 0.988. 

The best-fit LDGP model happens to be an LCDM 
model, using the current data, but it is more meaningful 
to give the 68% confidence limits on the mean values of 
the parameters. The 68% contour in Fig. shows that 
a broad range of non-LCDM models is consistent with 
the joint data constraints. In principle, we can calculate 
this range by marginalizing in turn over Q\ and Q m . In 
practice this is made difficult by the asymmetrical shape 
of the physical region in parameter space. The marginal- 
ized likelihoods are shown in Fig. 0] The f2 m likelihood 
has a near-gaussian shape, so that the median gives a 
good measure of the mean. Then we can calculate the 
limits by integrating to 34.1% of the area to the left and 
34.1% to the right. This leads to 



0.261 



i-0.019 
-0.018 ' 



(36) 



By contrast, the likelihood curve is strongly asymmet- 
rical, reflecting the fact that the unphysical region has a 
stronger effect on than on Q m . As a consequence, 



the mode and the median are significantly different, and 
we cannot extract a meaningful mean value for f^A- A 
conservative upper limit on Oa may be computed by in- 
tegrating from the mode until 68% of the area under the 
curve. This gives 



< 0.95 at 68.3% confidence . 



(37) 



If we use the Gold SN data [2(| instead of the SNLS 
data, then there are some interesting changes in the val- 
ues of the parameters, but the best fit lies still on the 
LCDM line. The Gold data favours a higher best-fit f2 m : 



Q m = 0.275 , il A = 0.725 . 



(38) 



X 



179.3 , x 2 P er degree of freedom = 1.142 .(39) 



For LCDM the x 2 P er degree of freedom is 1.135. The 
marginalized likelihoods for Q m and J1a are shown in 
Fig. We find that 

O m = 0.274±g:gil , (40) 
A < 0.87 at 68.3% confidence. (41) 

Thus the Gold data gives a stronger upper limit on Oa- 
Equations l|36|) . (|37|l and l|40|l. (|41|l summarize our 
quantitative conclusions about the observational con- 
straints on LDGP models. The 3 independent data 
sets allow for a significant range of these phantom-like 
braneworld models, with Q m values that are consistent 
with other observational tests, and with J7a values that 
imply significant screening of A. However, the LCDM 
model, with one less parameter, does provide a better fit 
to the data. The evidence for braneworld screening of A, 
via effective phantom behaviour, is therefore weak. This 
is consistent with results for general relativistic models. 
If one uses only supernova observations, especially the 
Gold data |20|. then there appears to be evidence for 
phantom dark energy. But it is misleading to ignore the 
CMB observations, and when these are imposed, the ev- 
idence for phantom behaviour is much weaker |2l| . 

It is important to stress that we have only imposed ob- 
servational tests on the background dynamics of LDGP. 
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FIG. 5: As in Fig. |1J but using the Gold supernova data. 



In regard to background dynamics, which arc determined 
by H(z), the LDGP model is indistinguishable from a 
GR model with dark energy whose equation of state is 
exactly w e g(z), as given by Eq. J2J. However, the GR 
model has no physical motivation, whereas the LDGP 
model gives a covariant and consistent physical meaning 
to the equation of state l|21l) . Nevertheless, one needs to 
break this degeneracy. The way to do this is to impose 
observational tests based on structure formation. 

Why have we not done this, i.e., computed the de- 
tailed CMB anisotropies and matter power spectrum for 
the LDGP model? The reason is that the density pertur- 
bation equations for LDGP are 5-dimensional, and their 
solution is still a formidable unsolved problem. An ap- 
proximate solution to the 5D equations in the DGP(+) 
models has been developed for small, sub-Hubble scales, 
by Koyama and Maartens |23j , confirming the intuitively 
motivated solution of Lue and Starkman [24|. This ap- 
proximate solution may also be applied to the LDGP case 
via the replacement r c — > — r c , as pointed out by Lue and 
Starkman |9(. But perturbations on scales near to and 
above the Hubble scale feel increasingly strong 5D grav- 
ity effects, and the approximation in Refs. p3.l2^| breaks 
down. The perturbation problem becomes strongly 5D 
and the Fourier modes satisfy partial rather ordinary 



differential equations. In the DGP(+) case, the CMB 
anisotropies and matter power spectrum have been pre- 
sented in Refs. . These papers describe valuable 
strategies for distinguishing DGP from GR dark energy. 
However, as explained by Koyama and Maartens |23|| . 
their quantitative results are unreliable, since they ignore 
5D effects on density perturbations (which turns out to 
violate the 4D Bianchi identity). 

Finally, we comment on the DGP ghost issue. In self- 
accelerating DGP(-I-) models, the asymptotic de Sitter 
state of the model suffers from a ghost j25j. However, in 
the LDGP case, which is based on the DGP(— ) branch, 
there is no ghost j2|J. 
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